Reconstructing the dark energy equation of state with varying alpha 
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The possibility of reconstructing the dark energy equation of state from variations in the fine 
structure constant is investigated for a class of models where the quintessence field is non-minimally 
coupled to the electromagnetic field. For given classes of couplings and quintessence interaction 
potentials, it is typically found that variations in alpha would need to be measured to within an 
accuracy of at least 5 x f -7 to obtain a reconstructed equation of state with less than a twenty per 
cent deviation from the true equation of state between redshifts and 3. In this case, it is argued that 
the sign of the first derivative of the equation of state can be uncovered from the reconstruction, thus 
providing unique information on how the universe developed into its present dark energy dominated 
phase independent of high redshift surveys. Such information would complement future observations 
anticipated from the Supernova Acceleration Probe. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Some recent observations of a number of quasar ab- 
sorption lines indicate that the fine structure constant, 
a = e 2 /He, was smaller than its present value by Aa/a — 
— 10~ 5 at redshifts in the range z ~ 1—3 PIQ- (See, how- 
ever, Ref. 0,0 for an independent analysis that does not 
support such a large variation in a.) Since this redshift 
range coincides with the epoch when the universe under- 
went a transition from matter domination to dark energy 
domination 0. . it is of interest to consider the possibil- 
ity that this change in the effective fine structure constant 
arises as a direct result of a non-trivial gauge coupling 
between the dark energy and the electromagnetic field 
strength 
(see also 

In this paper we consider classes of models where the 
dark energy in the universe is identified as a slowly vary- 
ing, self-interacting, neutral scalar "quintessence" field 

m us m m H3 ( - see also the reviews ) th&t is 

minimally coupled to Einstein gravity but non-minimally 
coupled to the electromagnetic field. The action is given 
by 




5' 



(i) 



where R is the Ricci curvature scalar of the metric 
g = detg^, k 2 = %-Km'p 2 and rap is the Planck mass. 
The Lagrangian density for the quintessence field is 

1 



-d^d^ -v(4>), 



(2) 



where V(4>) is the self-interaction potential of the scalar 
field, </>. The interaction term between the scalar field 
and the electromagnetic field is determined by 

C^ F = ~B F (<f,)F^F^ , (3) 



where is the electromagnetic field strength and 

Bf{4>) represents the gauge kinetic function that 
parametrizes the coupling between the scalar and vec- 
tor degrees of freedom. Cm represents the Lagrangian 
density for the ordinary matter fields and we assume this 
sector to be dominated by a barotropic pressureless fluid. 

The action Q is characterized in terms of two un- 
determined functions - the gauge kinetic function and 
the self-interaction potential. In general, these would 
be determined by the nature of the underlying parti- 
cle physics theory. For example, a generic exponential 
coupling of the type given in Eq. Q was introduced by 
Bekenstein pyi]. Exponential couplings between form- 
fields and scalar fields also arise generically in compact- 
ifications of string/M-theory to four dimensions, where 
the scalar field parametrizes the volume of the extra di- 
mensions. (See, e.g., Ref. 37] for a recent review on 
the theoretical motivation of varying fundamental con- 
stants.) 

The gauge kinetic function specifies the value of the ef- 
fective fine structure constant such that a = ao/ Bp((f>), 
where a subscript '0' denotes the present-day value. The 
potential of the field is related to the dark energy equa- 
tion of state, w<p = V4>l Pin where = <j) 2 /2 — V((j)) 

represents the pressure of the field, = <f> 2 /2 + V{<j>) is 
the energy density and a dot denotes differentiation with 
respect to cosmic time. 

A number of different approaches may be adopted 
when fitting models of the form (JIJ to the data. In prin- 
ciple, the functions {V((f>), Bp (</>)} would be determined 
within the context of a unified theory of the fundamen- 
tal interactions, such as string/M-theory. In this case, 
a direct approach would be to determine the region of 
parameter space consistent with observations once these 
two functions have been specified. This approach was 
effectively followed recently by Parkinson, Bassett and 
Barrow |l5|. who calculated the best fit parameters of a 
model for which the exact form of the gauge kinetic func- 
tion and the dark energy equation of state were assumed 
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a prion. 

We adopt an alternative approach in the present work 
by considering whether the quintessence potential or the 
gauge kinetic function can be reconstructed directly from 
observational data involving variations in the fine struc- 
ture constant and the dark energy equation of state. For 
a given w(z) and Aa/a = (a — ao)/ao, there always ex- 
ists a Bp((p) that would fit the data. Thus, a sufficiently 
accurate empirical determination of the equation of state, 
together with the evolution of a, would allow the gauge 
kinetic function to be reconstructed. 

On the other hand, if the gauge kinetic function 
alone is specified a priori (either through theoretical 
or phenomenological considerations), w(z) and Aa/a 
may no longer be viewed as independent variables, since 
they share a common origin through the rolling of the 
quintessence field. In effect, a consistency relation exists 
between these two quantities and an observational deter- 
mination of one would constrain the other. This implies 
that the study of the absorption lines in quasar spectra 
can in principle yield additional information on variations 
in the dark energy equation of state (and the correspond- 
ing quintessence potential). This is important, given that 
a determination of the redshift dependence of the dark 
energy equation of state directly from the luminosity dis- 
tance relations is difficult - the latter is determined by 
a double integral over the former and this can severely 
restrict the available information on the equation of state 
that can be extracted from observations [38|, |3!j ■ 

In this paper, we consider a linear dependence of the 
gauge kinetic function on the scalar field: 



B F (cf>) =l-C«fa-&) 



(4) 



where £ is a constant. This dependence may be viewed as 
arising from a Taylor expansion of a generic gauge kinetic 
function and is expected to be valid for a wide class of 
models when k(4> — 4>o) < 1 is satisfied over the range 
of redshifts relevant to observations, z « — 4. It then 
follows that the effective fine structure constant depends 
on the value of the quintessence field such that 



t depends 



Aa 

a 



a 



(5) 



Assuming that the mass of the scalar field effectively 
vanishes, tests of the equivalence principle imply that 
the parameter, £, is bounded by |C| < 10 -3 |9j. Bounds 
on variations in the fine structure constant arise from 
the Oklo natural nuclear reactor (|Aa/a| < 10~ 7 at red- 
shift z = 0.14 |4lLl42l | |56|). and the meteorite constraint 
(|Aa/a| < 1CT 6 at redshift z = 0.45 [i^V 

For consistency, neither the constraint arising from the 
Oklo natural nuclear reactor nor the meteorite constraint 
were considered in this work, although it is generally ex- 
pected that significant variations should be observed up 
to a redshift of order unity in quintessence models. How- 
ever, the model in Fig. [3| naturally satisfies the former 
bound at redshift z = 0.14. As discussed in Ref. liCJl. 



these bounds can be evaded through the existence of a 
form factor in the coupling £ with respect to the photon 
momentum. Such a form factor can result in changes in a 
at the level of atomic physics without leading to observ- 
able effects on nuclear phenomena. Moreover, it has been 
shown for a particular model in Ref. |25|. that if dark 
energy collapses along with dark matter this would nat- 
urally lead to a significant difference between the value 
of the fine structure constant in our galaxy and the one 
in the background. For these reasons, we decide to ex- 
plore also the models that do not satisfy the Oklo and 
meteorites bounds at low redshifts. 



II. RECONSTRUCTING THE EQUATION OF 
STATE: IN PRINCIPLE 

To proceed we consider a spatially homogeneous 
quintessence field propagating in the spatially flat 
Friedmann-Robertson- Walker (FRW) universe. We as- 
sume that the contribution of the electromagnetic degrees 
of freedom to the total energy density of the universe is 
negligible and consequently that the cosmic dynamics is 
determined by the scalar field and a background pres- 
sureless fluid (corresponding to dark and visible matter). 
It then follows that the cosmic dynamics is determined 
by the Einstein equation 



HH' 



and scalar field equation 



(pm + H 2 



subject to the Friedmann constraint 



H 2 = y(PM 



+ P<t>)- 



(6) 



(7) 



(8) 



where a prime denotes differentiation with respect to 
N = In a, H = a/a is the Hubble expansion parame- 
ter, the matter density is given by pu = Pa^Mo/o? and 
Po denotes the present value of the critical energy density. 
The term (dB p (</>) / d<j) F lw F ,tv ) containing the derivative 
of the gauge kinetic function was neglected in Eq. JJJ as 
its statistical average over the present Hubble radius is 
zero for photons 

Only two of the equations ©-JSJ are independent and 
the cosmic dynamics is fully determined once the func- 
tional form of the quintessence potential, V((f>), has been 
specified. The nature of the potential determines how 
the field evolves in time and the corresponding varia- 
tions in the fine structure constant are then determined 
by Eq. J3J). Moreover, the equation of state is defined in 
terms of the field's kinetic and potential energies: 



_ 4> 2 - 

W — - 

2 + 2V(<f>) 



(9) 
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and its dependence on redshift also follows given the form 
of the potential. 

The principle idea of the inversion procedure is that 
variations in the fine structure constant may be em- 
ployed within the context of this class of models to de- 
duce changes in the kinetic energy of the scalar field at 
higher redshifts. Eqs. ©-0 ma Y then be employed to 
determine the corresponding changes in the field's energy 
density, or equivalently, its potential, and hence the equa- 
tion of state from Eq. 0. This procedure is analogous 
to that employed in analysing the classical dynamics of 
a particle, x, moving in a one-dimensional potential well. 
The form of the well determines the particle's motion, 
and this can be represented as a trajectory in the phase 
space parametrized by {x,x}. Equivalently, the corre- 
sponding potential can be reconstructed once the appro- 
priate trajectory has been specified. 

We now develop the inversion procedure. Substituting 
Eq. (JHJ into Eq. (JJJ) yields a differential equation for the 
evolution of the energy density of the scalar field: 



a' = -( K 0')V + a~ 3 ) 



(10) 



where we have defined a = p,p/ po^Mo- The general solu- 
tion to Eq. IjlUII can be expressed in terms of quadratures 
with respect to the kinetic energy of the quintessence 
field: 

a(N) = e -I N dN(^'f x 



TV 



(11) 



where the integration constant ao is defined such that 
a(N = 0) = er = £10 o /^mo- 

The dark energy equation of state is given in general 
by w = — 1 — (m<7)'/3 and substitution of Eq. ifTUI) implies 
that 



w(N) 



-1 



IV2 



1 



1 



(12) 



Hence, the equation of state can be reconstructed once 
the redshift dependence of the first derivative of the field 
has been determined. It is important to emphasize that 
only the first derivative needs to be measured. If the 
gauge kinetic function is known, this dependence can be 
inferred directly from variations in the fine structure con- 
stant. 

It follows from the definition <j> 2 = p^ + that the 
equation of state l|12[) can also be expressed in the form 



3rL 



(13) 



where = k 2 p^ / (3H 2 ) . Differentiating Eq. (fT~5|> with 
respect to scale and substituting Eqs. ®), and © 
then implies that the first derivative (running) of the 



equation of state can be expressed directly in terms of 
first two derivatives of the scalar field: 



w' = 2(1 + w)-- + 3w 1 



(14) 



where ■ 



- (dw I dz)j a. An infinite hierarchy of expres- 



sions relating the n-th derivative of the equation of state 
to the (n + l)-th derivative of the field could be derived. 
Each represents a consistency relation between the equa- 
tion of state and variations in the fine structure constant 
once the gauge kinetic function has been specified. 

We now illustrate the above reconstruction procedure 
with a specific example where k(</> — 4>o) — AiV for some 
constant A. In this case, the integrals in Eq. (|llf> can be 
evaluated analytically: 



tf>0 



A 2 



V. 



Mq 



-X N 



A 2 



-3N 



(15) 



where A 2 = 3f^ (l + wq). The equation of state is then 
deduced by substituting Eq. (fT3f) into Eq. (JT32J: 

w(N) = (A 2 - 3) :•• - 



A ' exp ((A 2 - 3)JV)1 



w tt 4 

(16) 

Finally, the quintessence potential can be reconstructed 
by noting that V(N) = £ImoPo& — <t>' 2 H 2 /2 and employ- 
ing the Friedmann equation JSJ). We find that 



Be 



(17) 



where the mass scales A and B are positive-definite and 
given by 



A = - 



1 A 2 



2 3-A^ 



(18) 
(19) 



respectively. In the above example, it was assumed im- 
plicitly that the gauge kinetic function was such that the 
cosmological variation of Aa/a corresponded to a varia- 
tion in the scalar field of the form (f> oc N. The form of 
Bp(4>) was not specified. 



III. RECONSTRUCTING THE EQUATION OF 
STATE: IN PRACTICE 



In this section we consider the reconstruction of three 
dark energy equations of state by employing the method 
outlined in the previous section for a gauge kinetic func- 
tion given by Eq. (0J. The scalar field potentials have 
been investigated previously within the context of viable 
quintessence models [lJi Hji \Mm ■ The reconstructions 
are shown in Figs.^ CH These examples correspond to 
three different possible evolutions for w(z), namely, those 
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cases where it increases, decreases or oscillates with in- 
creasing redshift. The second and third examples are 
particularly important as they can not be reproduced 
with the parametrization employed in Ref. |l5j| . since in 
that work the equation of state is always increasing with 
increasing redshift. For the models studied in the present 
work, we have verified that k(4> — (f>o) < 1 over the ap- 
propriate range of redshifts, and this is consistent with 
the interpretation of Eq. Q as a lowest-order Taylor ex- 
pansion of a generic gauge kinetic function. Let us now 
describe the reconstruction process. 



A. Generating simulated data 

In testing the reconstruction procedure, it is necessary 
to first generate simulated data sets for the variations 
in a. This was achieved by specifying the functional 
form of the quintessence potential and numerically in- 
tegrating the field equations 10-© to determine the 
redshift dependence of both the equation of state and 
the quintessence field. The latter determines the corre- 
sponding variations in a from Eq. (0 once the coupling 
parameter, (, has been specified. We then generated 
the simulated data set, with associated error bars, for 
Aa/a based on the exact numerical solution. Specifi- 
cally, the data points are equally spaced in the redshift 
range z £ [0.2,4] at intervals of 0.2 and are normally 
distributed with mean £k(</> — 4>o). In each example, the 
value of £ was chosen so that the variations in — 4>o) 
resulted in changes in the fine structure constant of the 
order Aa/a ~ 10 -5 , as observed in the present QSO 
data 0. 



B. Fitting the data 

The reconstruction can then proceed by fitting the gen- 
erated data points to a polynomial function 



Aa 

a 



g!N + g 2 N 2 + ... , 



(20) 



where gi are constants. The result of these fits is shown 
in Fig. 0] Equations © and (|20|) map the variations in a 
as measured by QSO observations onto the correspond- 
ing variations in the scalar field for a given value of the 
coupling constant, (. The degree of the polynomial H2()(l 
employed in the fitting differs for the three cases because 
the underlying equations of state exhibit different levels 
of complexity. For the class of models we have studied 
we found that a polynomial of degree three provides gen- 
erally a good fit to the generated data. However, for 
models with an oscillating equation of state, one needs 
to increase the degree of the polynomial in order to ob- 
tain both a successful reconstruction and a reduced x 2 of 
order unity. 



C. Estimating £ 

The first and second derivatives of the field arc related 
to g such that <fi' = g' /C,K an d <j>" = g" /C K > respectively. 
In practice, therefore, the numerical value of the cou- 
pling parameter must be estimated empirically since the 
reconstruction via Eqs. I|ll[) and (|12|l requires the scale 
dependence of the quintessence field to be known. Sub- 
stituting Eq. (j20J into Eq. O implies that 



9 ' 



3 



(21) 



and it follows from Eq. that a numerical estimate for 
the coupling may be deduced given the present-day values 
of the quintessence field's energy density and the equa- 
tion of state wo, together with the variation, g , in the 
fine structure constant, as determined from QSO obser- 
vations. For example, given the typical values Q,p ~ 0.7, 
-w ~ 0.6 - 0.99, and g' Q ~ 10~ 7 - 10~ 5 , we find that 
C ~ 10~ 7 — 10~ 4 , in accordance with the values obtained 
in Ref. 01 > where specific quintessence models were stud- 
ied. This range of values of £ is compatible with bounds 
arising from tests of the equivalence principle which de- 
mand ici < io" 3 0. 

However, one source of uncertainty in the reconstruc- 
tion procedure is the uncertainty in the present-day value 
of the equation of state, wq. The latest measurements 
constrain this parameter within the range —1.38 < wq < 
—0.82 at the 95% confidence level, assuming a constant 
equation of state . This uncertainty generates an un- 
certainty in the value of the coupling £ and, in view of 
this, we allowed the equation of state to take a range of 
possible values, wq- More specifically, in Figs. and EJ 
the present value of the equation of state was chosen to 
be Wq = wq + (—0.1, 0, 0.1) where wq represents the cor- 
rect value as deduced from the numerical integration. In 
Fig- El on the other hand, the present value was chosen to 
be wq = w + (1 + wo)(0.9,0, —0.95), respectively, when 
moving downward in the figure. 

It follows from Eq. I|21|l that, when w « — 1, a small un- 
certainty in the value of wq can lead to a large uncertainty 
in the value of the coupling constant £ and consequently 
to distinct possible evolutions for the corresponding equa- 
tion of state. A typical case of w « — 1 arises when the 
scalar field undergoes oscillations about the minimum of 
its potential. On the other hand, an expression equiva- 
lent to Eq. (|14f> for the form of Bp ((f)) adopted in this 
paper is given by 



c 2 



w - 



2 1 g" 

3 Om g' 



(22) 



It follows, therefore, that in such cases a more accurate 
estimate of the magnitude of Q can be made if informa- 
tion on the present-day value of the first derivative of 
the equation of state is also available. We note that the 
quantity dw/dz is an observable parameter believed to 
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be within reach of future Snla observations from the Su- 
pernova Acceleration Probe (SNAP) |48|. 

We must note, however, that if £ was known on funda- 
mental particle physics grounds, the full reconstruction 
of the equation of state could be achieved without any 
need to normalize it to an independent result. 



D. Reconstruction results 

The results of the reconstructions are illustrated in 
Figs, n - El where Eq. 1)2 111 has been employed to esti- 
mate the value of the coupling £ for the different choices 
of wo by substituting w — > wq and g' — ► g\ . In each case, 
the dashed line in the figures represents the exact numer- 
ical solution of the equations of motion when Q^, = 0.7. 
The corresponding solid lines illustrate the reconstructed 
evolution for the different values of wq considered. An 
uncertainty in the present value of the equation of state 
of Swq w 0.1 is expected from the SNAP data. Hence, 
the lines with wq = wo ± 0.1 (in Figs.nand[3J) define the 
error band on the evolution of w arising from the uncer- 
tainty we will have on wq. It is worth emphasizing that 
SNAP will provide data within the range of redshifts be- 
tween and 1.7, whereas the QSO data can (in principle) 
provide us with information on the equation of state out 
to a redshift as high as z — 4. 

When investigating the sensitivity of the reconstruc- 
tion procedure on errors in the variations of a, we typi- 
cally found that in order to obtain a reconstruction, w{z), 
with no more than a 20% deviation from the true equa- 
tion of state, w(z), i.e., 



(for Wo = wq), one requires an observational determina- 
tion of Aa/a to within an accuracy of at least ~ 5 x 10 -7 
between redshifts and 3. Figure 0] shows how small the 
error bars on data points should be in order to obtain a 
reliable reconstruction at this order. This is an order of 
magnitude smaller than the expected sensitivity of the 
High Accuracy Radial Velocity Planet Searcher spectro- 
graph 8 {Aa/a) w 10" 6 . H^. 

IV. DISCUSSION 

Cosmological observations including high redshift sur- 
veys of type la supernovae and the anisotropy power 
spectrum of the Cosmic Microwave Background (CMB) 
indicate that the present-day value of the dark energy 
equation of state is bounded by —1.38 < wq < —0.82 at 
the 95% confidence level, assuming a constant equation 
of state |47| . Such bounds would be weakened for a wider 
class of models where the equation of state is allowed to 
vary, but at present there are only very weak observa- 
tional constraints on the "running" of the equation of 
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FIG. 1: Dashed line: evolution of the equation of state of 
the scalar field <j> with potential V{4>) = Vo exp((K0) 2 /2)/</> 11 
determined from numerical integration of the field equations 
|44) . Solid lines: evolution of the reconstructed equation of 
state for different possible values of the present-day equation 
of state (see the text for details). The function g{N) was 
fitted as a polynomial of degree 3 to a set of points normally 
distributed with a standard deviation of 1 x 10~ 7 . 

-0.7 ■ . ■ 1 



-0.8 
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FIG. 2: Dashed line: evolution of the equation of state of 
the scalar field <f> with potential V{4>) = Vb(exp(50K^>) + 
exp(0.8K<?!>)) determined from numerical integration of the 
field equations |4q|. Solid lines: evolution of the reconstructed 
equation of state for different possible values of the present- 
day equation of state (see the text for details). The function 
g{N) was fitted as a polynomial of degree 3 to a set of points 
normally distributed with a standard deviation of 5 x 10~ 7 . 

state, dw/dz — —aw' and it is not yet possible to 
distinguish such models from a cosmological constant. 
In this paper, we have investigated the possibility that 
further information on the redshift dependence of the 
equation of state can be deduced independently of high 
redshift surveys through observed variations in the fine 
structure constant. In principle, the reconstruction of 
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FIG. 3: Dashed line: evolution of the equation of state of the 
scalar field (j> with potential V = K- A e~ AK ' t ' [(k(/> - C) 2 + B)] 
determined from numerical integration of the field equations, 
where A = 8.5 and B = 0.93/yl 2 dill]. Solid lines: evolu- 
tion of the reconstructed equation of state for different possi- 
ble values of the present-day equation of state (see the text 
for details). The function g(N) was fitted as a polynomial of 
degree 5 to a set of points normally distributed with a stan- 
dard deviation of 5 x 10 -7 . 
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FIG. 4: Evolution of Aa/a for the three quintessence 
models studied in this work. The error bars represent the 
maximum uncertainty allowed in order to obtain a reli- 
able reconstruction of the equation of state with no more 
than 20% deviation from that deduced from the numer- 
ical integration. (a): V(<p) = Vp exp((K(j}) 2 /2)/0 11 0; 

(b): V(</>) = V o (exp(5OK0) + exp(O.8K0)) gg ; (c): V = 
k -a c -a^ _ c y. + B ^ ^ where A = g 5 and B = o.93/^2 

[Til E3| . The large error bar in the figures represent a typ- 
ical current uncertainty on Aa/a from QSO observations of 
6 (Aa/a) M 3.5 x 10" 6 . 



the equation of state is possible if the form of the gauge 
kinetic function that couples the scalar and electromag- 
netic fields is known. The advantage of a reconstruction 
of this type is that it yields information on the equation 
of state at redshifts significantly higher than the limited 
range accessible to SNAP (corresponding to z < 1.7). 

The primary question addressed in the present pa- 
per is how much information one could acquire on w(z) 
from variations in the fine structure constant alone. In a 
full reconstruction, one would employ all the data avail- 
able, from both supernovae surveys and measurements of 
Aa/a, and perform a full cross analysis between the dif- 
ferent data sets. However, one must also establish what 
can be learned from each data set independently. Indeed, 
this is a necessary and crucial step in the program we 
have outlined, precisely because w(z) and Aa/a are not 
independent as they share a common origin through the 
quintessence potential. As a result, information on varia- 
tions in the equation of state determined sep arately from 
supernova surveys (see, e.g. 0, |H |H 1|) and quasar 
surveys (as presented above) should be consistent. Es- 
tablishing an inconsistency would indicate that a reliable 
reconstruction could not be achieved and, furthermore, 
would immediately rule out this class of models, namely 
the form of Bp(4>), as a mechanism for correlating dark 
energy and variations in a. 

Figs.n~El m dicate that the reconstructions do yield in- 
formation on whether the running of the equation of state 
is positive or negative, at least out to a redshift z rj 3. Al- 
though the error in the normalization of wq implies that 
the uncertainties in the magnitude of the reconstructed 
first derivative may be large, the qualitative shape of the 
equation of state can be deduced, provided variations in 
a are determined to within an accuracy of 5 x 10~ 7 or 
better. We have performed equivalent analyses for the 
potentials considered in Section III over different regions 
of parameter space, as well as for other quintessence po- 
tentials, and have arrived at similar conclusions. 

Any information that can be extracted directly from 
observations on whether the equation of state increases 
or decreases with redshift is of importance. For example, 
if the equation of state increases with redshift (i.e. w 
moves away from —1), this implies that the field is slow- 
ing down as we approach the present day. On the other 
hand, the kinetic energy of the field is growing as the 
universe expands if w decreases with increasing redshift. 
This latter behaviour could correspond, for example, to 
a creeping quintessence scenario |33| , where the field has 
overshot the attractor value and has started to move only 
very recently. Thus, information on the first derivative of 
the equation of state provides us with unique insight into 
how the universe underwent the transition from matter 
domination to dark energy domination. 

For the general class of models defined by Eqs. (|TJ and 
(J5J , the qualitative behaviour of the equation of state can 
be deduced directly from Eq. (|21|l without the need to 
solve Eq. (|1L)|) if it is observed that g' 2 increases with 
redshift. Since Eq. (|21|l is valid over all scales and 
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is a decreasing function of redshift, it necessarily follows 
that the equation of state must have been larger in the 
past and this case would therefore rule out the possibility 
of a creeping quintessence scenario. On the other hand, 
for the case where g' 2 is a decreasing variable, we must 
proceed to solve Eq. QlOfl directly in order to gain further 
insight. 

Finally, we outline a complementary approach that 
may allow the equation of state and its derivatives to 
be deduced at a specific redshift. This approach corre- 
sponds to a perturbative reconstruction of the equation 
of state. It follows from Eqs. (|14|l and JSJ) that the first 
derivative of the equation of state at a given redshift 
can be directly determined if the corresponding values 
of {w,g',g"} are known. Higher derivatives can also be 
constrained if sufficient information on the corresponding 
derivatives of the fitting function g(N) is available. As- 
suming that the necessary constraints on the derivatives 
could be determined from QSO observations, the one re- 



maining free parameter would be the equation of state, 
or equivalently from Eq. , the density of the dark en- 
ergy. This parameter could in turn be deduced from the 
Friedmann equation JSJl if the Hubble parameter, H(z), 
were known and this could be found from the luminosity 
distance, cIl- 

= I(tti) <"> 

It would be interesting to explore this possibility further. 
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